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The Sn-induced y/Z x \/3 surface reconstruction on Si(lll) has been investigated by material- 
specific many-body calculations and by angle-resolved photoelectron spectroscopy (ARPES). This 
triangular surface system in the low adatom coverage regime is governed by rather localized dangling 
bond orbitals with enhanced electronic correlations and it is prone to exhibit magnetic frustration. 
We find a rather good overall agreement of the spectral function and its temperature-dependence 
between theory and experiment. Changes in the ARPES band topology in comparison to the 
density functional calculations can be explained as a spectral weight redistribution with respect to 
an additional symmetry which is not due to any geometrical change. This new symmetry corresponds 
to a magnetic order, which is found to be more complex than the canonical 120° anti-ferromagnetic 
order on a triangular lattice with nearest-neighbor coupling only. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.30.+h 



The adatom-covered surfaces of Si(lll) and Ge(lll) 
provide an excellent playground to study the competition 
and cooperation of geometrical frustration and electronic 
correlations in quasi-two dimensional triangular systems. 
The deposition of 1/3 of a monolayer of the group IV 
elements Pb and Sn renders these surfaces a/3 x -\/3 re- 
constructed [1-3]. These adatoms induce states in the 
semiconducting band gap, including surface bands which 
are well separated from the bulk bands. These states are 
dominated by highly localized dangling bond orbitals, 
and their narrow band nature leads to enhanced elec- 
tronic correlations in these systems. Moreover, long- 
range magnetic ordering is generally suppressed by mag- 
netic frustration inherently contained in such effectively 
triangular surface lattices [4]. 

In the a-phase of Pb/Ge(lll) [1] and Sn/Ge(lll) [2], 
a transition from a \/3 x ^/S to a 3 x 3 structure occurs 
upon cooling, which can be explained within a dynamical 
fluctuation model [5] . Pb or Sn atoms oscillate vertically 
at room temperature (RT), but their motions are frozen 
out at low temperature (LT), leaving two out of three 
Pb (Sn) atoms per unit cell at a different height than 
the third one. A similar structural transition might thus 
be expected to be observed on the surface of Sn/Si(lll). 
However, Sn/Si(lll) shows no 3 x 3 signature in low- 
energy electron diffraction (LEED) and scanning tun- 
neling microscopy (STM) experiments at both RT and 
LT, down to 6K [3, 6]. Careful analysis of photoelectron 
diffraction results further revealed that all Sn adatoms 
have the same bonding geometry [7] . 

In addition, in Sn/Si(lll) the observed valence-band 
photoemission spectra show a shadow band at P^g (see 
Fig. 2 for notations) and an approximate 3x3 peri- 
odicity of the overall surface band topology [3, 8]. Also, 
experimentally a clear conductance dip at the Fermi level 



was observed in scanning tunneling spectroscopy at low 
temperature [7]. The ground state, thus, is believed to 
be a narrow gap insulator. Calculations based on den- 
sity functional theory (DFT) in the local-spin density 
approximation (LSDA) failed to explain both facts, i.e. 
a metallic ground state without backfolded band features 
at Pyg is predicted [9, 10]. The insulating ground state 
can be explained by the strong electronic correlations fa- 
vored by the small bandwidth of the Sn/Si(lll) surface 
band. In LSDA -f U [4] and more sophisticated local 
density approximation (LDA) -I- many-body calculations 
based on the Hubbard model [11, 12], the ground state is 
confirmed to be insulating and identified as a narrow gap 
Mott insulator. However, the shadow band feature and 
the additional band symmetry observed in ARPES have 
not been explained consistently with respect to all the 
above mentioned experimental facts. In particular, the 
inherent spin frustrations contained in Sn/Si(lll) fur- 
ther elaborate the issue [12-14] with a remaining lack of 
knowledge on their interplay with electronic correlations. 

In this Letter, we provide a consistent explanation 
for the surface band topology and the shadow band to 
originate from electronic correlations rather than from 
structural aspects. We provide direct evidence from 
our temperature-dependent calculations that the approx- 
imate 3x3 symmetry observed in ARPES is indeed a 
consequence of magnetic correlations. However, the or- 
der associated with it is more complicated than the stan- 
dard 120 ° anti-ferromagnetic (AFM) spin-arrangement 
for triangular surfaces. 

Fig. 1 shows our experimental and theoretical situa- 
tion for the study of the Sn/Si(lll) system. Experimen- 
tally, clean surfaces were obtained by thermal desorption 
of the capping oxide at « 1250 °C from n-type Si sub- 
strates (p < 0.01 ficm) resulting in a sharp 7x7 pattern 



in LEED. Subsequently, 1/3 monolayer Sn was deposited 
on the substrates by electron beam evaporation. After 
an anneal at « 700 °C, the Sn/Si(lll)-vl5 x VS surface 
reconstruction was verified in high quality by LEED and 
STM, see Fig. 1(a). Theoretically, we construct a slab 
with six Si layers and saturate the Si bottom layer by 
hydrogen atoms. Sn atoms are placed on the top layer at 
equivalent T4 lattice sites, see Fig. 1(b). This is in agree- 
ment with the atomic structure derived from surface x- 
ray diffraction [15]. Therefore, we exclude any additional 
lattice superstructure, e.g., a 3 x 3 periodicity, in our cal- 
culations and expose our system to the electron-electron 
interaction U only. 
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FIG. 1. (a). STM image of the Sn/Si(lll) surface at T = 
300K (bias: -1.0 V; 0.5nA). (b). Top view of the DFT-LDA 
slab consisting of six Si-layers sandwiched between one top 
Sn-layer and one bottom H-layer (not shown in this figure). 
Sn atoms align in a plane. A 9-site Sn cluster is used in the 
DC A calculations (see text). 

First, we carry out ab initio DFT calculations based 
on LDA [12], which predict a half-filled metallic ground 
state of Sn/Si(lll), see Fig. 2(a) and (b). This is 
well in agreement with preceding LDA calculations [9- 
11]. Next, we project the Sn-related surface band onto 
the maximally localized Wannier basis and construct a 
single-band Hubbard model. The dynamical cluster ap- 
proximation (DCA) [16] is used with the continuous-time 
quantum Monte Carlo method [17] to solve this model in 
the paramagnetic phase. We consider a 9-site cluster in 
our calculations, as shown in Fig. 1(b). The surface 
Brillouin zone (SBZ) is divided into 9 equal-area sectors. 
In each of them the electron self-energy becomes a con- 
stant, i.e., i;(fc, iw„) ^ 'E{K,iujn). Our LDA -I- DCA 
calculations reported here are the first to examine the 
spectral weight redistribution in momentum space in the 
thermodynamic limit for Sn/Si(lll). In what follows, 
we address the theoretical photoemission spectra, which 
are related in the usual way to the imaginary part of the 
single-particle Green's function. The latter is calculated 
following a recipe contained in Ref. [18]. 

Fig. 2 contains the experimental and theoretical pho- 
toemission spectra as a function of position in the SBZ. 
ARPES experiments were carried out at the electronic 
structure factory (ESF) endstation of beamhne 7.0.1 at 
the Advanced Light Source (ALS), which provides sam- 
ple cooling down to 10 K and is equipped with a 6-axis 
goniometer. Photoelectrons were detected with a Sci- 
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FIG. 2. (a), (b) LDA-derived band structure along the two 
major high symmetry directions FK^ and TM ^ in com- 
parison with (c), (d) corresponding photoemission spectra 
(T ~ 60 K, = 130 eV) and (e), (f) the spectral function 
as obtained from LDA -f- DCA calculations. In (c) and (d), 
the contour lines are shown as a guide to the eye. White ar- 
rows in (c - f) indicate additional band maxima, which are 
not present in the LDA-derived band at the same energies. 
The four yellow circles around in (c - f) highlight the 
shadow bands observed in both ARPES and in the LDA -|- 
DCA calculations, which are also missing in the LDA-derived 
band. 



enta R4000 spherical analyzer with energy resolution set 
to 25 meV throughout all measurements. In the theoret- 
ical spectra calculations, we set U ~ 0.66 eV, which is 
slightly above Uc ~ 0.60 eV for the metal- insulator tran- 
sition (MIT) in this system [12]. In a related system, i.e. 
Sn/Ge(lll), earlier studies reported two distinct surface 
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FIG. 3. (a). LDA band structure (solid lines), backfolded LDA bands (dashed lines) with respect to 3 x 3 SBZ boundaries, 
and ARPES band map (false color plot) with intensity maximum (red dash-dotted line) indicated as guide to the eye. In order 
to achieve a better agreement with the experimental band edge at K we have to scale the LDA-band, which is shown as a 
solid line in Fig. 3(a), by a factor of 1.3. This scaling effect has essentially been captured by the LDA + DCA calculations in 
Fig. 2. This factor should not be interpreted as the effective mass of this surface band, since Fermi liquid theory fails to explain 
the shadow band. (b). The 1st \/3 x \/3 SBZ of the Sn/Si(lll) surface, which is divided into three different sectors (A, B, 
C) in the DCA approximation, see text for more details. The black and violet hexagons correspond to the SBZ of a 3 x 3 and 
2\/3 X 2\/3 superstructure, (c). Average particle density in each momentum sector as a function of interaction, reflecting the 
weight transfer from sector Kc to Ka- 



states in ARPES [5, 19]. Sn atoms were believed to stay 
at two different adsorption sites [20-22] in this system. 
However, in the Sn/Si(lll) surface band, such a spUt- 
ting is very small at low temperature and even absent at 
room temperature [3, 8]. We do not want to study such 
small energy differences in our theoretical spectral func- 
tion, which contains the uncertainty from the analytical 
continuation. However, it should be noted that quantum 
Monte Carlo calculations have shown that strong elec- 
tron repulsions can further split the lower Hubbard band 
[23, 24]. Thus, whether or not the two peak structure in 
that system should be attributed to a geometrical change 
deserves a more careful study. Moreover, it was shown in 
a recent theoretical work that electron correlations alone 
can induce structural transformations in elemental iron 
[25]. Here, we will not address such additional effects. 
For Sn/Si(lll) we compare ARPES and theory by fo- 
cusing on the band topology and its k-dependent spectral 
weight. 

Two ARPES k-space line scans along TK ^ and TM ^ 
are shown next to the corresponding LDA -f DCA results 
in Fig. 2(c) to (f). Rather good overall agreement can 
be detected in this comparison. Along both TK ^ and 
FM^, we observe a shadow band in ARPES around 
F^, as indicated by yellow circles in Fig. 2(c) and (d). 
The evolution of this band clearly shows an additional 
hand maximum at a position close to along TK ^ 
and in vicinity of if 3 along TM ^ indicated by the white 
arrows. M3 and are high symmetry points of a 3 x 3 
SBZ. The appearance of the additional band maximum, 
which is absent in the LDA results in Fig. 2(a) and (b), 
modifies the spectral dispersion from ^/i x •\/3 to an ap- 
proximate 3x3 symmetry. In the LDA -I- DCA calcula- 
tions shown in Fig. 2(e) and (f), both the shadow hand 
and its spectral evolution as a function of momentum are 



well reproduced. The shadow band is clearly visible in 
theory and slightly more pronounced than in the experi- 
mental spectra. The overall agreement between ARPES 
and the LDA -I- DCA calculations, and especially the ap- 
pearance of this shadow band, represent strong evidence 
of many-body effects in this system. In contrast, a struc- 
tural origin is rather unlikely since in our calculations all 
Sn atoms are located at equivalent lattice sites within the 
same atomic layer. Thus, the approximate 3x3 symme- 
try cannot result from any structural change, which is 
essentially in accordance with the already mentioned ab- 
sence of a surface-band splitting in the Sn/Si(lll) system 
at low temperatures. 

Profeta and Tosatti suggested that the 120°-AFM or- 
dering and consequent folding might be the origin of the 
shadow band [4]. In what follows, we want to demon- 
strate that it is very likely that this system is magneti- 
cally short range ordered. However, a strict 120°-AFM 
order cannot fully explain the surface band topology, es- 
pecially the energy range of the shadow band and the 
position of the band maximum which is not located ex- 
actly at K-i along the Fi\f ^3 direction. According to 
Profeta and Tosatti, we backfold the original LDA-band 
into the 1st SBZ with respect to the 3x3 SBZ bound- 
ary, which then corresponds to a 120°-AFM order. In 
Fig. 3(a), the original LDA-band is plotted as a solid 
line along the VKMY ^ directions, the dashed lines are 
the folded bands. The high symmetry points of the 3x3 
SBZ are labeled as T^^M^^K^. Three different colors 
on the LDA and the backfolded bands are used to indi- 
cate different momentum sectors, as those in Fig. 3(b). 
These three inequivalent momentum sectors are derived 
from the symmetry preserved in the 9-site cluster DCA 
calculations, which are labeled as sector A, B and C. The 
high symmetry points, i.e., Fy^, M ^ and K are con- 
tained in sector A, B and C, respectively. The red dashed 
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line represents the Fermi surface (FS) of the LDA band, 
which is completely contained in sector B. Evidently, in 
the 3x3 SBZ a band is located around at i? < 0, 
which is obviously back-folded from K^, see Fig. 3(a). 

First of all, our calculations strongly support the exis- 
tence of a magnetic order and the consequent band back- 
folding. In Fig. 3(c), the average particle numbers n{k) 
at each momentum sector are shown as a function of the 
interaction U. n{k) was directly calculated in the LDA 
+ DCA. It relates to the spectral function A{k^ E) by 
n{k) = J^^dEA{k,E)f{E,T), where E denotes the en- 
ergy and f{E,T) the Fermi function. We found qualita- 
tively a different behavior of n{k), when k e Ka,Kb or 
Kc- n(Kj^) monotonically grows with the increase of the 
interaction, while n{Kc) behaves exactly the opposite. 
In contrast, h^Kb) stays almost constant while varying 
U. In sector B, with the increase of U, the quasiparticle 
peak at the Fermi level gradually looses its weight until 
a charge gap opens. The constant value of n{KB), there- 
fore, strongly indicates that the total spectral weight in 
sector B does not change, however, the spectral weight 
lost at the quasiparticle peak transfers to the lower and 
upper Hubbard bands within this sector. The constant 
value of u^Kb) is a strong indication of the Mott type of 
the MIT. Moreover, an increasing U results in a spectral 
weight transfer from sector C to sector A. In the U = 
limit, there is no intensity at Fy^ for energies E < 0, 
giving the almost zero value of n{KA)- For higher values 
of U, part of the spectral weight around transfers 
to F^, resulting in the increasing/decreasing behaviors 
of n{KA)/n{Kc)- Thus, what we observe from Fig. 3(b) 
mainly reflects the spectral weight transfer from to 
F^, which strongly supports the band back-folding pic- 
ture. 

However, a strict 3x3 symmetry cannot fully ex- 
plain the shadow band we observed. A close compari- 
son of ARPES and the folded LDA bands reveals that 
the shadow band at F^ is not at the same energy as 
the band at as it would be if the magnetic or- 

der was 120°-AFM. The shadow band stays at higher 
energies than that at K^. Thus, the magnetic SBZ 
of Sn/Si(lll)--y3 X \/3_R30° can only be approximate 
to 3 X 3. It reflects that the magnetic order derived 
from our ARPES and calculations is actually close to 
but different from the classic 120°-AFM for triangular 
systems. This is partially due to additional hopping pro- 
cesses inherently contained in Sn/Si(lll) as compared to 
the ideal triangular model with nearest-neighbor hopping 
only. In our previous study, we found that the inclusion 
of the next-nearest-neighbor hopping in triangular lat- 
tices changes the spin-susceptibility peak-position from 
near to in the 1st SBZ which indicates a mag- 
netic order change from 120°-AFM to a row-wise (RW) 
type AFM [12]. This is equivalent to the magnetic or- 
dering in the Mn/Cu(lll) surface, which can be effec- 



tively described by a triangular Heisenberg model with 
higher order exchange interactions [26]. The superposi- 
tion of three equivalent spin arrangements of RW-AFM 
order, which possess a 2^/3 x 2a/3 magnetic SBZ , can 
further lower the total energy and is thus favored. M3 
of the 3x3 SBZ is also a high symmetry point of the 
2V3 X 2^/3 SBZ. However, as can be seen from Fig. 3 
(b), -^2^3 of 2a/3 X 2-^3 is close to but different from K3 
of the 3x3 SBZ, which essentially explains the agreement 
of the band maximum position with M3 along TM ^ and 
the discrepancy with along TK^ . The current surface- 
band topology study and the spin susceptibility calcula- 
tions in our previous work [12] coincide with each other, 
and both point at the magnetic order to be of RW-AFM 
type. To this end, a spin symmetry-broken many-body 
calculation and a spin-resolved STM study are highly de- 
sirable for further understanding of this adatom system. 
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FIG. 4. Comparison of the angle-integrated ARPES spectra 
(color filled area) with the local density of states calculated 
from the LDA -I- DCA at (7 = 0.66 eV (solid black lines) for 
different temperatures. 

In Fig. 4, we show a comparison of the temperature 
dependent angle-integrated photoemission spectra with 
the local density of states (LDOS) calculated from the 
LDA + DCA. The experimental curves are obtained by 
integrating the ARPES results along the F^g — -^^3 ~^ 
K ^ — Fyg directions. Experimentally, we observe a 
MIT between T = 100 K and 150 K. The angle-integrated 
spectral weight at the Fermi level becomes nonzero with 
increasing temperature. This is in agreement with the 
many-body description of a Mott MIT. On the other 
hand, higher temperature has almost no effect on the 
spectra far away from Ep. We find a marginal change on 
the spectra for energies below — 0.6eV. In contrast, for 
increasing temperatures spectral weight is transferred to- 
wards the Fermi level, driving this surface system from 
insulator to metal. Theoretically, for [7 — 0.66 eV we ob- 
tain an overall good agreement with the experiments for 
different temperatures. At all temperatures represented 
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in Fig. 4, the part of the LDOS close to the Fermi level 
coincides well with its experimental counterpart. 

In summary, we have shown that the key features of 
a triangular adatom system, realized by the Sn/Si(lll)- 
V3 X \/3i?30° surface, can be qualitatively explained by 
strong electronic correlations. By assuming a planar con- 
figuration of the Sn atoms, we find a good overall agree- 
ment between experiment and the LDA + DCA calcu- 
lations. A temperature dependent MIT is found, which 
closely coincides with the Mott description for this sur- 
face system. We find strong evidence for a spectral weight 
transfer from the momentum region around to F^, 
indicating the existence of a magnetic order in this sys- 
tem. The additional symmetry observed in ARPES can 
then be understood as a band back-folding with respect 
to a new magnetic ordering. This should stimulate fur- 
ther studies on the magnetic properties of this and related 
systems. 
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